Pseudospin polarized quantum transport in monolayer graphene 
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Monolayer graphene with an energy gap presents a pseudospin symmetry broken ferromagnet with 
a perpendicular pseudomagnetization whose direction is switched by altering the type of doping 
between n and p. We demonstrate an electrical current switching effect in pseudospin version of 
a spin valve in which two pseudoferromagnetic regions are contacted through a normal graphene 
region. The proposed structure exhibits a pseudomagnetoresistance, defined as the relative difference 
of resistances of parallel and antiparallel alignments of the pseudomagnetizations, which can be 
tuned to unity. This perfect pseudomagnetic switching is found to show a strong robustness with 
respect to increasing of the contact length, the effect which we explain in terms of an unusually long 
range penetration of an equilibrium pseudospin polarization into the normal region by proximity 
to a pseudoferromagnet. Our results reveals the potential of gapped graphene for realization of 
pseudospin-based nanoelectronics . 

PACS numbers: 73.23.-b, 72.25.-b, 72.80.Vp, 85.75.-d 



I. INTRODUCTION 

Graphene, the two dimensional layer of the carbon 
atoms with honeycomb lattice structure, has attracted 
a great deal of attention as a new promising material for 
nanoelectronics, since its experimental realization a few 
years agcP^^l, Most of the peculiar properties of graphene 
is the result of its massless Dirac spectrum of the low- 
lying electron-hole excitations, which in addition to the 
regular spin appear to come endowed with the two quan- 
tum degrees of freedom, the so called pseudospin and 
valley. The pseudospin represents the sublattice degree 
of freedom of the graphene's honeycomb structure, and 
the valley defines the corresponding degree of freedom in 
the reciprocal lattic^^IlHll, The effect of these additional 
quantum numbers has already been exploited by anoma- 
lous features of several quantum transportphenomena in 
graphene, including quantum Hall effeclpnS, conductance 
quantizatiorfiSlj Klein tunnelin^SHU! and quantum shot 
noise'™'. 

Interestingly, the pseudospin and valley degrees of free- 
dom in graphene have been proposed separately to be 
used for controlling the electronic devices in the same 
way as the electron spin is used in spintronic and quan- 
tum computing. Rycerz et al^^^^ demonstrated an elec- 
trostatically controlled valley filter effect in graphene 
nanoribbons with zigzag edge which can be used for re- 
alizing valley valve structures in valleytronics (valley- 
based electronics) applications. On the other hand a 
pseudospin-based version of a spin valve has been pro- 
posed in bilayer graphene, where the pseudospin is deter- 
mined by the relative amplitude of the wave function on 
the two layer^2ni_ A bilayer pseudospin valve consists of 
two connected neighboring regions whose pseudospin po- 
larizations can be tuned by application of gate voltage 

In this paper, we study the possibility of realizing pseu- 
dospintronics in monolayer graphene within the scatter- 
ing formalism. The possibility of an interaction driven 
spontaneous breaking of the pseudospin symmetry, which 
can lead to the realization of pseudomagnetic states 



in monolayer and bilayer graphene, has been studied 
recentljI^H, Here, we demonstrate that the monolayer 
graphene with a gap in its electronic spectrum and an 
appropriate doping presents a pseudospin symmetry bro- 
ken ferromagnet, with a finite pseudospin magnetization 
oriented vertically to the graphene plane. The magnitude 
of the pseudomagnetization (PM) depends on the chem- 
ical potential and its direction is switched by changing 
the type of doping (electron n or hole p). 

Based on the above observation, we propose a 
nonmagnetic pseudospin valve which consists of two 
pseudoferromagnetic (PF) regions separated by a non- 
pseudomagnetizaed normal (N) graphene of length L 
(shown schematically in Fig. [ijj). The PM direction in 
each region can be tuned independently by means of elec- 
trical gates, which allows for switching between parallel 
and antiparallel configurations. We find that the pro- 
posed PF/N/PF spin valve exhibits a relative difference 
of the electrical resistance in the two parallel and an- 
tiparallel states of the pseudomagnetizations, which can 
be remarkably large in analogy to the giant magnetore- 
sistance (GMR) in magnetic multilayer J^. When the 
chemical potential of the system is tuned to be close to 
the energy gap (/i ~ A), the pseudospin valve effect can 
be perfect PMR takes the value 1 for appropriate choices 
of the length L. More importantly, we show that the 
perfect pseudospin valve effect can be reached even in 
higher chemical potentials 3> A by applying an appro- 
priate bias voltage. We further demonstrate the unusual 
proximity effect at pseudoferromagnet-normal junctions 
(PF/N, PF/N/PF). We find that an equilibrium pseu- 
dospin polarization is induced into the N region with a di- 
rection which is precessing around the axis normal to the 
PF/N interface, and an amplitude which decays slowly 
in the distances of few Fermi wave length Xp from the 
interface. This is in clear contrast to the induced magne- 
tization in ordinary ferromagnet-normal metal junctions, 
which decays exponentially within Xj^. 

This paper is organized as follows. In Sec. |ll) we intro- 
duce pseudoferromagnets (PF) and use them to study the 
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FIG. 1: (Color online) (a) Vertical pseudomagnetization per 
electron PM^/N of the gapped graphene layer versus chemical 
potential p ( is scaled to the energy gap A), (b) Schematic 
illustration of the proposed pseudospin valve in monolayer 
graphene: The left and right regions are pseudoferromagnets 
(PF) and the intermediate region is a normal graphene (N) 
without a band gap. (c-d) Profile of pseudomagnetization 
vector PM inside the two PFs (blue) and the N region of 
length L = Xf (pink) for two configurations of (c) parallel 
and (d) anti-parallel, when ~ A. 



pseudospin valve effect in monolayer graphene PF/N/PF 
junction. Sec. |III| is devoted to the investigation of the 
proximity effect in graphene PF/N junctions. Finally, we 
present the conclusion in Sec. llV] 



II. PSEUDOSPIN VALVE 

The pseudospin valve structure we study is shown 
schematically in Fig. [TJa. It consists of two PF re- 
gions with tunable direction of PM which are con- 
nected through a normal (non-pseudomagnetized) layer 
of length L. Such a structure can be realized in a 
graphene sheet with a sizable gap in its electronic band 
structure. There are several methods to open an en- 



ergy gap in the band structure of graphene. A sce- 
nario is placing graphene on top of an appropriate sub- 
strate which breaks the graphene sublattice symmetry 
and generates a Dirac mass for charge carriers. The band 
gap opening is ob served in epitaxially grown graphene 
on a SiC substrat^^HHI and a hexagonal boron nitride 
crystaP^, and also on metallic surfaces with a boron ni- 
tride buffer layei-^. The energy band gap engineering can 
be also achieved through doping the gra phen e with sev- 
eral molecules such as CrOs, NH3, H20^^^. A gapped 
graphene with a chemical potential ^ close to its energy 
gap A acquires a PM of sizable magnitude and an ori- 
entation which switches with changing the sign of the 
chemical potential (see Fig. [I^) . A highly doped gapped 
graphene has a vanishingly small PM and behaves as a 
non-pseudomagnetized (normal) graphene. 

To study quantum transport in the pseudospin valve 
structure within the scattering formalism, we first con- 
struct the quasiparticle wave functions that participate 
in the scattering processes. We adopt Dirac equation of 
the form 



Hib ^ (e + 



(1) 



where in the presence of an energy gap H = vp{(T.p) + 
A(Tz is the two-dimensional Dirac Hamiltonian with 
p — ~ih^ the momentum operator in the x-y plane 
{vp — 10^ m/s represents the Fermi velocity) and cr = 
{cxi o-yjO'z) the vector of the Pauli matrices operat ing in 
the space of two sublattices of the honeycomb lattic^^^^. 
The two-dimensional spinor has the form ip = {'iPa,4'b)i 
where the two components give the amplitude of the wave 
function on the two sublattices and e is the quasiparticle 
energy. 

For a uniform gapped graphene region the solutions of 
the Dirac equation Eq. ([T]) are two states of the form 



e± _ ±ik^x iqy 



e 



(2) 



for conduction band electrons of n-doped graphene and 



e± _ „TiKx iqy 



(3) 



for valance-band electrons of p-doped graphene, 
at a given energy e and transverse wave vec- 
tor q with the energy-momentum relation 
ec,v ^ ±[-M + a/A2 + {hv\kc,y\y]. Here 
ac,v = aTCsm{hvq/ ^/{e±~f^j^~^^A^) is the angle 
of propagation of electron which has longitudinal 
wave vector fcc,D = {^vf)~^ y/{e ± ^)2 — A^ cos a and 
4>c,v — arcsm/i(A/-^/(F±~/ip^-~A2). The two propaga- 
tion directions of electron along the x axis are denoted 
by ± in . 

The pseudospin of such states for conduction (valance) 
band electrons of n- (p-)doped graphene is obtained as 
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FIG. 2: (Color online) Pseudomagnetoresistance (PMR) of 
the pseudospin valve versus the length of the N region {L/Xf) 
for different values of n/A. 



where k± is the unit vector normal to the electronic wave 
vector fcc,t)- 

As can be seen from the above equation, the existence 
of a band gap makes the pseudospin to have a component 
perpendicular to the electronic momentum in the plane of 
the graphene sheet. The in-plane and out of plane com- 
ponents of the pseudospin depend on (e -|- fj.) /A, which 
can be tuned to unity to make the pseudospin vector to be 
oriented perpendicular to the sheet. Increasing (e-t-/i)/A 
leads to the decrease of the out of plane component such 
that it goes to zero when e + fi ^ A. 

The total PM of the gapped graphene is calculated 
by summing the expression Q over all the wave vectors 
k = {k,q), 



(5) 



from which we find that PM only has an out of plane 
component which depends on (e + ^)/A. Fig. [TJi shows 
the behavior of the out of plane component of PM per 
electron PMz/N as a function of /i/A at zero tempera- 
ture (T ^ 0). It is seen that for a* - A, PMJN takes 
its maximum value PM^/N = 1 , while increasing /.t/A 
leads to the decrease of PM^ such that it goes to zero 
for highly doped gapped graphene (/i ^ A). For a given 
value of the orientation of PM changes by changing 
the sign of /J./ A. This implies that the PM vectors of 
gapped graphene regions with different type (n or p) of 
dopings are oriented antiparallel versus each other. 

We note that the pseudospin polarization of a gapped 
graphene corresponds to a difference in the electronic 
charge densities of the two triangular sublattices, which 
in turn produces an in-plane electrical polarization of the 
density Pn,p — "^kPn-p- electron with a wave vector 
fc contributes a dipole moment which is calculated from 



±(1 



') e aAB, 
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where a^s is the lattice vector oriented from A to B site. 
We see that the electric dipole moments are oriented from 
A to B (B to A) sites for n- (p-)doped gapped graphene. 
This correspondence between PM vector and the density 
of the in-plane electrical polarization can be used for an 
experimental measuring of PM. 

With the above found behavior of PM, the realization 
of the pseudospin valve of Fig. [TJd seems quite feasible. 
The configuration of the pseudospin valve can be changed 
from parallel to antiparallel by fixing the type of doping 
of one region and changing the type of the doping in 
the other region. The size of the pseudospin valve effect 
is determined by the extent in which the conduction of 
the anti-parallel configuration is suppressed (similar to 
the spin valve effect). The pseudomagnetoresistance of a 
pseudospin valve is defined as 



PMR- 



Rap — Rp Gp — Gai 
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where Rp(AP) is the resistance of the parallel (anti- 
parallel) configuration and Gp(^AP) is the corresponding 
conductances, which can be calculated from the Lan- 
dauer formulaP^ 



(8) 



Gp{AP)=9o / \tp(AP)\ COS a da, 



where Tp^AP) — I^p(ap)P is the transmission probabil- 
ity of electrons through the pseudospin valve in parallel 
(anti-parallel) configuration and = Ae^ /h is the quan- 
tum of conductance. 

We have calculated the amplitudes of the transmission 
tp(AP) by matching the wave functions of three regions of 
the left PF, N region and the right PF (signed by 1,2, and 
3, respectively) at the two interfaces, a: = and x — L. 
The solutions of Eq. ([T]) in the three regions for parallel 
(P) and antiparallel (AP) configurations are as follow 
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-0^1^^^ are the wave functions of Dirac equation for in- 
coming and outgoing electrons of n- (p-)doped graphene 
sheet with (without) a gap; r is the reflection coefficient 
in the left PF and a and b are the coefficients of the 



states 'ip. 



'e± 



Matching these solutions at the interfaces. 



we obtain the transmission amplitudes as 

~ C + D + e'"3-'#'3 {B + (e2»fe2i _ l)e-''"2) ' 
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where A = 4cosaiCOSQ!2 ea;p(i[2iisa. + q;2 + ^2-^]), 
B = exp{iai + icjyi) [exp{2ia2) + exp{2ik2L)], C ~ 
[1 — exp(2ifc2L)] ea;p(i[Q;i + (12] + 0i) and _D = 1 + 
ea;p(2iQ!2 + 2ik2L). 

Finally from the calculated expressions of Gp and Gap 
via Eq. ^ and using Eq. ([t]), we obtain PMR as a 
function of the chemical potential and the length of the N 
region. Fig. [2] shows dependence of PMR on L/\f{\f ~ 
hvp/^) at T = and zero bias voltage V ~ and for 
different values of n/A. We have taken /zi = ^2 = ImsI = 
/u. We observe that the pseudospin valve effect can be 
perfect {PMR = 1) for /i ~ A. For these values of /i, 
PMR shows an oscillatory behavior with L/Xp, with an 
amplitude which takes the value 1 for some ranges of the 
length. We note that this perfect pseudospin valve effect 
of the monolayer graphene is more robust with respect to 
an increase of the length of the N region, as compared to 
the similar effect in a bilayer graphene pseudospin valve 
structurc'^'^. 

The amplitude of PMR decreases by increasing /u/A 
and tends to the constant value of PMR = 1/3 for a 
highly doped structure with A. This residual PMR 
is the difference in the resistance of a n-p graphene struc- 
ture with that of a uniformly (p or n) doped graphene 
with the same which is present even in the limit 
PM — > of a non-pseudomagnetized structure. 

We have found that the perfect pseudospin valve ef- 
fect can be resumed by applying a bias voltage V of an 
appropriate amplitude to the valve. This is shown in 
Fig. [3]i,b. In Fig. |3}i, PMR is plotted versus L/\p for 
fi/A = 2 and different eV/A. As it is seen, by increasing 
eV from the amplitude of PMR increases and the oscil- 
latory behaviour is suppressed, such that PMR reaches 
the perfect constant value 1 when eV — A. Fig. [Sja 
shows plot of PMR as a function of eV/A for different 
values of L/Xp, which shows that the pseudospin valve 
becomes perfect as eV — 1, independent of the value 
of L/Xp. Applying bias voltage is such that there is no 
band crossing for p-doped excitations {eV < fi~A). This 
is similar to the case of having PF regions with different 
chemical potentials, where = eV ± fj, and /x„ > \fip\. 



III. PROXIMITY EFFECT IN PF/N 
JUNCTIONS 



The above found strong robustness of the pseudospin 
valve effect with respect to increasing the length of N 
contact should arise from a strong pseudomagnetic cou- 
pling between the two graphene PF regions. This itself 
can be due to a long range penetration of pseudospin po- 
larization into the N region by proximity to PF regions. 
To analyze this in detail, we study proximity effect in 
PF/N, PF/N/PF junctions. We start with a single PF/N 
junction in a graphene sheet in the x-y plane, where the 
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FIG. 3: (Color online) (a) Plot of PMR versus L/Xp for three 
values of bias voltage eV/A = 0, 0.5, 1 (b) and the behavior of 
PMR versus eV/A for three values of L/Xp, when ^ = 2 A. 



region a; < (PF region) has a uniform PAI oriented ver- 
tically to the sheet and the region x > (N region) is in 
the normal state. We calculate PM in PF and N regions 
using Eq. ([s]) , and by considering the contribution of the 
pseudospin of all incident electrons from left and right 
regions that are scattered from the junction. The total 
pseudomagnetization vector PM will be the sum of PM 
vectors contributed from electrons incident from the left 
and the right regions: 



pm: 



PM, 1 PM 

— -i -\ —} 

Nj 2^ TVj NJ ^' 

PMf)=^<^(fc) >v,^,,^„, 



(12) 



where j denotes the PF(N) region. 

The resulting profile of PM across the PF/N junction 
is demonstrated in Fig. |4] for /i ~ A. As it is seen, a 
nonzero PM is induced in N region (A = 0) which ro- 
tates around the normal to the junction (x axis) with x. 
The perpendicular component PM^ oscillates as a func- 
tion of X with a period of order A^-, and shows only a 
weak decay in the scale of Af- While the in-plane com- 
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FIG. 4: (Color online) PM of the PF/N junction when /i ~ A: 
The profile of PM in PF (blue) and N regions (pink) (a) and 
the position dependence of the PM components (b). The 
inset of Fig. |4j3 shows the magnitude of PM. 



ponents PM^^y vanish inside the PF regions, they are 
produced at the PF/N interface and are penetrated into 
the N region. This is originated from the difference be- 
tween the pseudospin of electrons coming from PF side 
(Eq. ([4]) ) and that of electrons coming from N side (Eq. 
^ with A = 0), which results in a nonvanishing aver- 
aged over momentum directions pseudospin polarization. 
The y component PMy shows an oscillatory behaviour 
with X similar to PM^. Interestingly, PM^ is uniform 
inside N, which considering the decay of the other two 
components, implies that PM at the points in N far 
from the junction is uniform and oriented perpendicular 
(along X axis) to the PM in the connected PF. This 
unusual proximity effect can be explained in terms of re- 
flectionless Klein transmission of electrons which incident 
normally to PF/N interfacd^l^J'^. We note to the unusu- 
ally long range penetration of the proximity induced PM 
inside the N region, which is in contrast to the ferromag- 
nct /normal metal junction (FN), in which the induced 
magnetization decays over short interatomic distances. 

The above analysis of the proximity effect in PF/N 
junction can be extended to the pseudospin valve geome- 
try of Fig. [TJa. The profile of PM orientation in different 
regions of the PF/N/PF junction is indicated in Fig. [l] 
for parallel (c) and antiparallel (d) cases when L = Xp 
and /I ~ A. PM is perpendicular to the x-axis and un- 
dergoes rotation across the N contact in a way that in P 
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FIG. 5: (Color online) Plot of PM dependance on the length 
of the N region L/\f for PF/N/PF junctions, where fi ~ A. 



and AP cases PMy and PM^, respectively, shows change 
of signs at the middle of N region (a: = L/2). 

Furthermore, we obtain that the magnitude of PM 
inside the N region is constant with x for both of parallel 
and anti-parallel configurations. Dependence of the mag- 
nitude of PM on the length of the N region L and for 
/i/A ~ 1, is shown in Fig jsjfor both of parallel and anti- 
parallel configurations. We see that they have a periodic 
behavior with L. As is expected, in the limit L — > pseu- 
domagnetization of the parallel configuration goes to the 
constant value of a PF layer \PM\/N — 1, while it tends 
to zero for anti-parallel configuration (n-p junction). 



IV. CONCLUSION 

In conclusion, we have investigated realization of 
pseudospin polarized quantum transport in monolayer 
graphene with an energy gap in its Dirac spectrum. We 
have demonstrated that a gapped graphene is a pseu- 
dospin symmetry-broken state which exhibits a pseu- 
domagnetization PM oriented normal to plane of the 
graphene. The magnitude of PM depends on the ra- 
tio of the chemical potential to the energy gap /i/A 
and its direction is switched by changing the type of 
doping between n and p. Based on this observation, 
we have proposed a pseudospin valve in which a non- 
pseudomagnetized normal region connects two pseudofer- 
romagnetic (PF) regions, whose magnetization alignment 
can be controlled by altering the type of their doping. 
The suggested pseudospin valve exhibits a pseudomagne- 
toresistance PMR, defined as the relative difference of the 
resistances in parallel and antiparallel alignments, which 
for ~ A can be tuned to unity by appropriately ad- 
justing the contact length L. We have shown that this 
perfect pseudomagnetic valve effect with PMR = 1 is 
preserved even for very large lengths L ^ Xp- Although 
PMR decreases by increasing /i/ A in the absence of a bias 
voltage, a perfect switching at large chemical potentials 
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can be resumed by applying an appropriate bias voltage. 

In order to explain the robustness of the pseudospin 
valve effect with respect to increasing of the contact 
length, we have further studied the proximity effect in 
PF/N junctions. We have found that an equilibrium 
pseudospin polarization can be induced in the normal 
graphene over a very large length Xp. The induced 
pseudomagnetization vector PM undergoes a damped 
spatial precession around the normal to the PF/N jimc- 
tion and tends to be uniform along the normal at the 



large distances Xp from the junction. 
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